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A FEASIBILITY STUDY OF DETERMINING PHYSICAL PROPERTIES OF 
THE MARTIAN ATMOSPHERE BY USE OF SOLAR OCCULTATION 
AS SEEN FROM AN ARTIFICIAL SATELLITE 

By H. Andrew Wallio and James R. Williams 
Langley Research Center 

SUMMARY 

A feasibility study to determine the physical properties of the Martian atmosphere 
by the occultation of sunlight as seen by an artifical satellite has been conducted. Five 
model atmospheres have been used in this study. The effects of refraction, Rayleigh 
scattering, and limb darkening and finite sun size have been included in the mathematical 
model for visible light passing through the atmosphere. The study has shown that the 
technique of using the occultation of the sun as seen by an artificial satellite does provide 
a means of determining the physical properties of the Martian atmosphere provided that a 
differential correction technique to match observational data with theoretical data is devel- 
oped. This report also shows the importance of using the shorter wavelengths of visible 
light since the observed effects are most pronounced in this region of the spectrum. 

INTRODUCTION 

Currently, there are several unmanned vehicles planned for the exploration of the 
planet Mars. These vehicles will be used to gather data concerning the physical proper- 
ties of the planet. One of the primary areas of research will be the determination of the 
constituents of the Martian atmosphere and their distribution. This report presents a 
relatively simple technique which should be considered as a means of determining such 
atmospheric properties as surface pressure and scale height. The technique employed 
utilizes the observed light -intensity variation from the sun as the light passes at grazing 
incidence through the planetary atmosphere. 

The theory describing the observed optical effects that occur during an eclipse of a 
satellite by a planet with an atmosphere is given by Link. (See ref. 1.) The application 
of occultation to the observed light intensity of the moon of a planet is discussed by Rakos 
(ref. 2) and Weisberg (ref. 3). The effects of scattering (refs. 4 and 5), absorption 
(ref. 6), and the finite size of the sun (refs. 1 and 3) must also be taken into account when 
the lower regions of the atmosphere are viewed. 



Since this paper is primarily a feasibility study, only the primary factors deter- 
mining the variation of light intensity have been included. These factors are refraction 
through the atmosphere, limb darkening, finite sun size, and Rayleigh scattering. Light- 
intensity-variation curves are presented for each of these effects considered both singly 
and collectively. The curves presented are for five model atmospheres and for a range 
of four wavelengths of visible light. Also a short discussion of the effect of the distance 
of the satellite from the planet is included. 

This technique utilizing the occupation of sunlight as seen by an artificial satellite 
can be used as an independent data type to determine scale height and surface pressure or 
to provide a check on any other independent data source. The advantage of occupation 
data is that by assuming an attitude control system and cruise orientation similar to that 
used on the lunar orbiter and Mariner spacecraft, one spacecraft axis would nominally be 
pointing toward the sun, and sensors properly mounted would require no additional atti- 
tude control. Also the use of an orbiter results in a number of observations and would 
allow a study of the effects due to clouds or dust since there should be a large number of 
observations that would be essentially cloud and dust free. 

SYMBOLS 


B 


normalized coefficient in limb darkening 


D 


distance from trajectory to planet (see fig. 2) 


ds, ds’ differential spacing of adjacent light rays (see fig. 3) 

F(D,x) integrated intensity due to finite size of sun 

f(D,x) intensity of light at point x a distance D from planet due to a point source 

at center of sun 


f 3 


g 

h 

I 

L(y) 


coefficient of depolarization of jth gas 
acceleration due to gravity, 3.7 m/sec^ 


R 0 T 0 

pressure scale height, , km 


mg 


intensity of light 


limb darkening of sun 


2 


path length, km 


factor determining solar shape and limb darkening 
average mass of gas, kg/mole 
index of refraction 

atmospheric pressure, bars (l bar = 10^ N/m^) 

gas constant per unit mass 

universal gas constant 

radius of planet, km 

temperature, °K 

variable limit of integration 

coordinate of satellite trajectory, km 

integration variable over solar disk 

altitude above surface of planet, km 

half-angle subtended by sun as seen from satellite, radians 
Rayleigh scattering coefficient, km - * 
temperature gradient, km"^- 

angle between line of sight to sun and sun planet line, radians 
first variation 
refr activity 
Loschmidt's number 


angle between minimum altitude and radius to satellite, radians 



X 


wavelength of light, km 


p. constant of integration 

v fraction by volume of ith gas 

p atmospheric density, kg/km^ 

r optical depth 

co deviation of light ray from straight line, radians 

Subscripts: 

e earth surface conditions 

i incident or initial conditions 

j summation index over number of gas constituents 

0 Martian surface conditions 

1 values at minimum altitude 

ASSUMPTIONS AND ANALYSIS 

The construction of a mathematical model in order to perform a feasibility study for 
the use of solar occultation data as a method of determining the physical properties of the 
Martian atmosphere has required certain assumptions which are outlined as follows: 

(1) The sun, planet, and satellite are all coplanar (fig. 1) 

(2) The planet is a sphere of radius r 

(3) The distance into the geometric shadow x is small compared with the distance 

from the observer to the center of the planet (See fig. 2.) 

(4) The radius of the planet is much greater than any characteristic height in the 

atmosphere 

(5) The atmosphere is a perfect gas of uniform chemical composition in hydrostatic 

equilibrium 
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(6) The refractivity is small compared with unity 

(7) The change in temperature is small over an altitude distance of one scale height 

The temperature profiles assumed are a least-squares linear fit, between 0 and 100 km 
above the Martian surface, to the temperature curves of the Martian atmosphere engi- 
neering models used by the Viking Project. The basic properties of the atmosphere are 
given in table I. These models as presented in table I provide the basis for this feasi- 
bility study and are the only models discussed herein. It is further assumed that the 
source is a finite disk and sufficiently far from the planet so that the incident light may 
be considered to be plane waves. The final assumption is that multiple scattering and 
Mie scattering are negligible. The possible importance of Mie scattering is discussed in 
reference 7. 

In this paper the three major effects upon the variation of light intensity are refrac- 
tion through the Martian atmosphere, limb darkening of the sun, and Rayleigh scattering 
due to the gas molecules in the atmosphere. 

The first effect under consideration is that of differential refraction which causes 
the divergence of two infinitesimally separated parallel rays as they traverse an atmo- 
sphere whose index of refraction is a function of the altitude. Figure 3 shows two inci- 
dent parallel rays of separation ds which traverse the atmosphere and emerge with a 
separation ds' which is greater than ds. Thus an incident beam with an intensity 1^ 
and an infinitesimal separation ds has a flux proportional to 1^ ds. Since the flux 
entering the tube shown in figure 3 must be equal to the flux emerging at the observation 
point, it follows that 


I = 



( 1 ) 


Therefore the intensity at the point of observation will necessarily be less than the inci- 
dent intensity. The spreading of the beam and the derivation showing the dependence of 
the spreading on the constituents and physical properties of the atmosphere are given in 
appendix A. 

Another effect which causes a diminution in the observed light intensity is due to the 
limb darkening of the sun. The light which reaches the satellite is a sum of contributions 
from the various parts of the disk of the sun, each passing through a different layer of the 
Martian atmosphere. The effect is a spreading of the light curve and a smoothing of the 
discontinuities due to the fine structure of the atmosphere. Implicit in limb-darkening 
effects is the consideration of the fact that the sun has a finite size. These two effects 
are included in this occultation analysis, and their derivations are included in appendix B. 
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The third major effect taken into account is that of Rayleigh scattering which is the 
scattering of light from the gas molecules in the atmosphere. Rayleigh discovered that 
for scattering in which the wavelength of the incident electromagnetic radiation is long in 
comparison with the dimensions of the scatterer, the scattering is dependent upon the 
fourth power of the wavelength. Radiation traversing the atmosphere suffers loss due to 
absorption and scattering. For this paper the effects of absorption were considered to be 
negligible in comparison with Rayleigh scattering. The intensity of a parallel beam of 
light decreases, because of the extinction process, by a factor of 

I = Ije"^ (2) 

where /3 is a local extinction coefficient derived by Rayleigh and l is the path length. 
The complete derivation of the dependencies of the physical properties of the atmosphere 
on /3 and l is given in appendix C. 

RESULTS AND DISCUSSION 

In order to facilitate the discussion, it is convenient to classify the various model 
atmospheres under study. Table I presents a list of the five models, used in this feasi- 
bility study, along with their constant physical properties. Throughout the remainder of 
this text, these models are referred to according to the numbers given in table I. 

Since this study is a feasibility study, a circular orbit of 8700 km was chosen arbi- 
trarily for the radius of the satellite orbit. For the effects which are dependent upon the 
distance from the planet, a separate discussion is included but, in general, the 8700-km 
circular orbit is of primary importance. 

In order to present a deeper insight into the effect of refraction, Rayleigh scattering, 
and limb darkening, each of these effects has been studied independently and in combina- 
tions with each other. 

The effect of refraction only is the first of these areas to be studied. Since for a 
low-density gas, the refraction is essentially independent of wavelength, wavelength is not 
included in figure 4. Figure 4 is a plot of the ratio of light-intensity variation as a func- 
tion of the distance into the geometrical shadow of the planet (see fig. 2) for the effect of 
refraction only. These curves are for a satellite in an 8700-km circular orbit. The 
significance of scale height and, in this case, surface temperature is illustrated in the 
difference in light -intensity variation in figures 4(a) and 4(b) and in the difference in 
light -intensity variation in figures 4(d) and 4(e). The only differences in physical proper- 
ties of figures 4(a) and 4(b) are surface -temperature conditions in which there is a 100° K 
difference and a small difference in the slope of the temperature profiles, likewise for 
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figures 4(d) and 4(e). Figure 4(c) represents the light -intensity variation for the so- 
called "most probable model." Note that in figure 4 all the light -intensity curves are 
terminated at approximately the same point. From equation (A29), it can be shown that 
for D = 8700 km, the point at which the light intensity drops to zero is approximately 
26.6 km. As would be expected, model IV gives the greatest variation in light-intensity 
distribution since it has the lowest surface temperature and the highest surface pressure 
of any of the models. Note in connection with figure 4 that the light intensity drops very 
slowly until it reaches a point at which it disappears completely. This effect is precisely 
the expected effect from considerations of refraction only from a finite source. 

The effect of limb darkening and finite size of the sun as discussed in appendix B 
is independent of wavelength. In order to best illustrate the effect of limb darkening and 
finite sun size on the light intensity, a computation was made for a planet with no atmo- 
sphere. Figure 5 illustrates this effect for a satellite in a 8700 -km circular orbit. As 
can be seen from figure 5, the light -intensity ratio drops from 1.0 to zero gradually over 
a distance of approximately 56 km. Since the effects of considering the finite size and 
the limb darkening of the sun are dependent upon the angle which the sun subtends as 
seen from the satellite, figure 6 illustrates this distance-dependent effect. As can be 
seen from figure 6, the greater the satellite orbital radius, the slower the light-intensity- 
ratio curves drop to zero. 

The next area of interest in this study is the effect of Rayleigh scattering. Figure 7 
illustrates the effect of Rayleigh scattering which is proportional to the inverse fourth 
power of wavelength. In order to illustrate the effect of Rayleigh scattering better, fig- 
ure 7 includes only Rayleigh scattering effects. In figure 7 it can be seen that the shorter 
wavelengths for a given model tend to suffer the greatest reduction in intensity. A note 
concerning figure 7 is that at approximately 26.6 km, the light is physically blocked from 
the satellite by the planet. Therefore, at this point the ratio of light intensity drops to 
zero. The effect of scale height can readily be seen in the difference between the curves 
in figures 7(a) and 7(b) where the major difference in models is the 100° K difference in 
surface temperature, and the small difference in the slope of the temperature profiles, 
likewise for figures 7(d) and 7(e). Figure 7(c) again presents the Rayleigh scattering 
effect for the most probable model. The effect of surface pressure as well as scale 
height is apparent in the difference between figures 7(a) and 7(e), and in the difference 
between figures 7(b) and 7(d). Note in comparing these figures that the vertical scale is 
not the same for all models. The effects shown in figures 7(a) and 7(e) and in figures 7(b) 
and 7(d) reflect the dependence of the Rayleigh scattering coefficient on density as can 
also be seen in the derivation in appendix C. 

Now that the magnitude and the importance of the single effects considered in this 
report have been discussed, it is necessary to combine the effects in order to obtain the 
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results of this study. For the sake of brevity, the summation of the effect of refraction, 
Rayleigh scattering, and limb darkening and finite sun size is hereafter referred to as the 
combined effects. In order to obtain the combined effects at any given point x, one 
simply takes the value of the light-intensity ratio at that point from each of the curves 
illustrating the single effect (figs. 4, 5, and 7) and forms their product. This product 
produces the resultant light-intensity variation for figures 8 to 11. 

Figure 8 is a plot of the variation of the light-intensity ratio with the distance of the 
satellite into the geometric shadow for a received light intensity of 6500 angstroms with 
the combined effects. In order to facilitate the discussion, it is convenient at this point 
to use the scale height h which is a surface-pressure scale height, as a factor for dis- 
cussing each of the models. The scale height associated with each of the models is pre- 
sented in table I. Figures 8(a) and 8(b) are models which have the same physical proper- 
ties except for a difference of 100° K in surface temperature and a difference in the slope 
of the temperature profile. This difference is reflected in a scale-height difference of 
approximately 5 km. For figures 8(a) and 8(b), the surface pressure of which is low 
(4 mb), the effect of a scale-height difference is small and only visually detectable in the 
region between -40 and -20 km into the shadow. Basically, the same observation may be 
made for figures 8(d) and 8(e); the only difference is a scale-height difference of approxi- 
mately 6.8 km. However, for figures 8(d) and 8(e), the surface pressure is 20 mb, and 
the difference in the light -intensity variation in the -40 to -20 km region of the shadow is 
more perceptible. Figure 8(c) is the light-intensity variation for the atmospheric model 
which reference 1 calls the most probable model. The differences among figures 8(a), 
8(c), and 8(d) are easily perceptible even though their scale heights are within 2.6 km of 
each other. This difference is attributed chiefly to their surface-pressure differential. 

Figures 9 and 10 present the same combined effects as were presented in figure 8 
except that figure 9 is for wavelengths of 5500 angstroms and figure 10 is for a wavelength 
of 4700 angstroms. One of the basic reasons for including figures 9 and 10 is for a refer- 
ence in case a cross plot involving wavelength is desired by the reader. 

Figure 11 is a plot of the variation of light-intensity ratio with the depth into the 
geometrical shadow for a wavelength of 3500 angstroms. Figures 11(a) and 11(b) as were 
figures 8(a) and 8(b), 9(a) and 9(b), and 10(a) and 10(b) are presented for models the only 
difference of which is in scale height. However, it should be noted that in figures 11(a) 
and 11(b), the difference between these two figures in the -40 to -20 km range is more 
perceptible. This effect is very obvious for figures 11(d) and 11(e). Again the scale 
height differs by 6.8 km; however, the difference in light curves is much more predomi- 
nant at the shorter wavelengths. This condition is a result of the Rayleigh scattering 
being inversely proportional to the fourth power of wavelength. 


8 



The results of this study can best be discussed with the aid of figures which contain 
the superposition of some previously presented curves. These curves illustrate the dif- 
ference in the light -intensity variation for the extremes in visible light and illustrate the 
ease with which model atmospheres can be distinguished when the lower wavelengths of 
visible light are used. 

Figure 12(a) is a superposition of figures 8(a), 8(c), and 8(d) which are for a wave- 
length of 6500 angstroms and a surface pressure of 4 mb, 9 mb, and 20 mb, respectively. 
The scale-height differential between models I and DI is 1.2 km with a surface-pressure 
differential of 5 mb. From figure 12(a), it can be seen that this value produces approxi- 
mately a 2-percent maximum difference in light curves. On the other hand, the scale- 
height difference between models in and IV is 1.4 km with a pressure differential of 
11 mb. This value produces a maximum variation in the light-intensity curve of approxi- 
mately 6 percent. Thus, if the scale height were known, the light sensor would require 
0.5-percent accuracy to determine surface pressures within 2 to 3 mb at 6500 angstroms. 

Figure 12(b) is a superposition of figures 8(d) and 8(e) showing that at 6500 ang- 
stroms, a scale-height differential of 6.8 km results in a 3-percent maximum difference 
in light-intensity variation at a constant surface pressure of 20 mb. 

Figure 13(a) is a superposition of figures 11(a), 11(c), and 11(d) for a wavelength of 
3500 angstroms and a surface pressure of 4 mb, 9 mb, and 20 mb, respectively. With a 
surface -pressure differential of 5 mb and a scale -height differential of 1.2 km between 
models I and HI, it can be seen from figure 13(a) that the light -intensity variation differs 
by a maximum of 11 percent. For a scale-height difference of 1.4 km and a surface- 
pressure differential of 11 mb for models in and IV, figure 13(a) shows a 33-percent 
maximum difference in light -intensity variation. Furthermore, there is a maximum dif- 
ference of 44 percent in the light-intensity variation for models I and IV as shown in fig- 
ure 13(a) for which the scale-height differential is 2.6 km and the surface-pressure dif- 
ferential is 16 mb. 

Figure 13(b) is a plot of figures 11(d) and 11(e) for the wavelength of 3500 angstroms 
and a surface pressure of 20 mb. The difference in scale height is 6.8 km which leads to 
a 11 -percent difference in light-intensity curves. Thus, the advantage of performing this 
experiment with the shorter (3500 angstroms) wavelengths is adequately illustrated in fig- 
ures 12 and 13. 

There are some general observations that can be made concerning figures 8 to 13. 

It appears that the surface pressure can be approximated by following the slope of the 
curves in the -60 to -20 km range, the higher rate of change of the slope in the region of 
the knee of the curve denoting the lower pressures. For curves with the same surface 
pressure and with the same wavelength of received light but different scale heights, the 
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curve with the lower scale height intersects the curve with higher scale height from above 
as shown in figure 13(b). These facts are mentioned because they would provide a basis 
for the formulation of a differential correction method to determine scale height and sur- 
face pressure from observational data. On the other hand, if one is provided with scale 
height by some external measurement, this technique will provide a method of determining 
surface pressure. 


CONCLUDING REMARKS 

This study has included the effects of refraction, limb darkening, finite sun size, 
and Rayleigh scattering. These effects have been considered singly and collectively. 
Singly, both refraction and Rayleigh scattering show effects of scale-height differences; 
however, limb darkening and finite sun sizes are the dominant effects in the light -intensity 
variation, and the predominance of scattering over refraction effects on the light curves 
has been shown. In using the occultation of the sun as seen by an artificial satellite as 
means of determining the physical properties of the Martian atmosphere, it has been 
shown that the most predominant effect, a 44 -percent variation in light intensity, is 
observed when visible light of short wavelengths (3500 angstroms) is used. 

This technique of occultation can provide the means of determining the physical 
properties, scale height, and surface pressure. However, the technique would depend 
upon the development of a differential correction technique to match the observed curves 
with the theoretical curves. However, this study has shown that the technique is a fea- 
sible one for determining the physical properties of the Martian atmosphere and warrants 
further study and consideration. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 17, 1969. 
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APPENDIX A 


ATMOSPHERIC REFRACTION OF LIGHT 

The vertical distribution of pressure in a fluid can be expressed in terms of the 
hydrostatic equation 


dp = -gp dZ (Al) 

Since an atmosphere consists of a compressible gas, its density is a function of altitude 

P = p(Z) (A2) 

Gravity too is a function of altitude above the planet’s surface; however, for this paper it 
was assumed to be a constant. 

In order to establish the variation of density with altitude, it is assumed that the 
atmosphere obeys the ideal gas law: 


p = pRT 


where 



Temperature is assumed to vary linearly with altitude over the pertinent range 


(A3) 


(A4) 


t = t 0 (i + rz) 


(A5) 


where T 0 is the surface temperature and T 0 T is the slope. Substituting equation (A5) 
into equation (A3) yields 


P 


P _ P 
RT RT 0 (1 + TZ) 


(A6) 


Substituting equation (A6) into equation (Al) yields 


f P dp = r z dZ 

Jp o P J 0 RT 0 (1 + TZ) 


(A7) 


11 



APPENDIX A 


The integral in equation (A7) gives 


P = P Q exp 




where 


1 

i + rz 


* 1 - 


rz 


(A8) 


has been used. 

The index of refraction for a rare gas may be written as 


n(Z) = 1 + e(Z) 


(A9) 


where e(Z) is the refractivity as given by the Lorenz -Lorentz formula which states that 
the refractivity is proportional to the density of the gas. Therefore, 


e(Z) =l2p(Z) 


Substituting equation (A6) into equation (A10) yields 


e(Z) = e 0 (l - TZ)exp 



(A10) 


(All) 


where e 0 may be evaluated from standard temperature and pressure conditions at the 
earth’s surface by using the relation 


e 


o 



(A12) 


The path of a ray of light passing through a planetary atmosphere can be found from 
Fermat’s principle 


6 



ds = 0 


(A13) 


where 


ds 2 = dZ 2 + (r + Z) 2 d6> 2 

The Euler -Lagrange equation which yields the minimum for equation (A13) is 
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n(Z)(r + Z ) 2 jj| 


1 


1 + (r - Z) 2 (f' 2 


1/2 


= Constant = / u. 


dZ 


At the point of closest approach, Zj, — = 0, therefore 

M = (l + e j) (r + Zj) 

Substituting equation (A15) into equation (A14) and integrating yields 


e < z > - £ 


i 1 + e l) 

( r + z l) 

1 dZ 

' + Z) 

[(1 + e) 2 (r + Z) : 

2 - 1 

(l + ei ) 2 (r + Zj] 

l 2 ~ 

1/2 


From figure 14, it can be seen from the geometry that 

U> . 7 T 




and 


r 

2 Jr 


(r + Zj) dZ 


Z 1 (r + Z)[(r + Z) 2 - (r + Z]) 2] 1 / 2 


Then it follows that 


“( z o = 2 £ 


r + Zi , 

(2rZ - 2rZ j + Z 2 - Zj 2 ) 1 

72 ^ 

1 

| 1 + + e £ 1 ) (r + z)2 -( r + Zl ) 

-|l/2 

2 

r + Z ' 

(2rZ - 2rZj + Z 2 - Zj 5 

1)1/2 

-J 


■ >dZ 


If it is assumed that when the integrand is large, 

(a) e j «- 1 

(b) r » |z - Zj 

(c) r » Z 
e - e x 


(d) r 


Z - Zi 


« 1 


(A14) 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 
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ll■■lllllllllllllllll 


then, 


f °° l 

? 1 + e -‘i r y 1/2 _ i 

- z ^ 1 / 2 

l z - Zl ) 


dZ 


e /. 3 6 1 " € 


By using approximation (d), equation (A20) becomes 

W ( Z l) a ]/l I Zi ( z 7^372^ + fz^Z! 

By using equation (All) in equation (A21), 


r + . . 


. dZ 


00 




1 - exp 


/Z - Zi 


z - z 


hi 


exp 


+ r- 


(z - Zi) 3 / 2 (z - Zi) 1 / 2 


i j (1 - exp 

h l / 3 

+ 4 e l r 




2 \ 


(z - Zl ) 


5/2 


/ 


where 


h i = 


RTj 


Integrating equation (A22) yields 

a, ( z i) = € i/^ 1 + ~ y l + (^' 


Generalizing equation (A23) for any altitude allows the subscript 1 to be dropped. 

The ratio of incident to received intensity is equal to the ratio of ds' to ds or 

I _ ds 
Ii ds' 

From figure 15, it can be seen that approximately 

I ds 


Ij ds + D do> 


14 


(A20) 


(A21) 


(A22) 


(A23) 


(A24) 


(A25) 


(A26) 


I 



APPENDIX A 



(A27) 


(A28) 


(A29) 


which is the result used in equation (1) of the text. For further discussion see refer- 
ences 3 and 4. An alternate approach to the same form is given in reference 8. 
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APPENDIX B 


LIMB DARKENING AND THE FINITE SIZE OF THE SUN 


The light from the sun which reaches a satellite is the sum of contributions from 
various parts of the sun, each ray passing through a different layer of the atmosphere. 
This condition results in a spreading of the light curve and a smoothing of the disconti- 
nuities due to the fine structure in the atmosphere. 

To account for this effect, assume that the distance D in figure 16 is a constant 
and the distance x is a variable function of the direction of the light ray emanating from 
the sun. 


If f(D,x) is the intensity at point x for a point source located at the center of the 
sun, the integrated intensity due to the finite size of the sun is of the form 


F(D,x) = f f(D, x+aDy)Z(y) dy (Bl) 


which is given in reference 2. The half -angle a subtended by the sun as seen from x 
and l(y) is a function which takes into account both the circular shape and the limb 
darkening of the sun. 

By assuming that the brightness of a segment of the sun is given by the product 
l(y) dy where dy is the segment width, the total light intensity of the sun is 



(B2) 


Let L(cos y) denote the factor giving the intensity of solar limb darkening where y is 
the angle between the line of sight to the sun at x and the normal to the sun; then in 
terms of x and y 


cos y 




(B3) 


where x and y as they appear here are normalized to the radius of the sun. The 
function l( y) can be expressed as the fraction of the light intensity received at dx due 
to the solar element dy. This condition results in 



APPENDIX B 


where the term in the denominator is a normalization factor as given in reference 2. By 
following reference 3 and assuming that L(cos y) is of the form 

L(cos y) = B + (1 - B) cos y (B5) 

equation (B4) becomes 

Uy) = 2 W X - y 2 + - B )( x - y 2 ) m 

This expression is substituted into the final form for the light intensity in appendix C. 


17 



APPENDIX C 


RAYLEIGH SCATTERING 


Radiation transversing a planetary atmosphere suffers losses in intensity due to 
absorption and to scattering in other directions. This effect is called extinction. Since 
the Martian atmosphere is rare, absorption has been neglected in this report. 

Atmospheric scattering will reduce the intensity of the incident light as given by 

I = Ije" T (Cl) 


where r is the optical thickness of the atmosphere and is given by 



(C2) 


/3 is the Rayleigh scattering coefficient, and d l is the differential path length. From 
reference 5 j3 is given by 


But 



(C3) 


(C4) 
(C 5) 


-2Z 

[ nj (Z) - l] 2 « (e oj ) 2 e h ° (C6) 

From the geometry (fig. 14) 

l 2 ~ (r + Z) 2 - (r + Zi) 

and in the region where ^ p dZ is large, r » Z + Zj so that l 2 ~ 2r(z - Zj.) and 
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cxp (© =exp E^ z v) 


(C7) 


Since, the number density t) q approaches 
becomes 


V o eXp l'2^ 



the optical depth 


T 


-Zl 


-7.2 




Collecting constant terms yields 


(C8) 



The constant terms correspond to the scattering coefficient at the surface, therefore, 
the optical depth reduces to 


r = 


p 0 fiirrh 0 expl- 



(Cll) 


where 

The primary effect of Rayleigh scattering is to vary the initial slope of the intensity 
curve as a function of the wavelength. Thus there should be a distinct intensity curve for 
each wavelength measured from the satellite. By combining the equations for refraction, 
Rayleigh scattering, and finite solar size, the final form for the equation of the light 
intensity as a function of x is 




exp l-j3 0 J2vrh 0 exp (- \ v 

^ = ° ^ j l( y) dy 


*)’ 


(° < z < z max) 


(Z = 0) 


(C12) 
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where 


and 


x + aDy = -Z + w, 


° D ( 1+ t 


rh„ - |rz1e h " 


(Cl 3) 




( is ^M i+ f n, °)' x ]) 
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TABLE L- SUMMARY OF CONDITIONS AT MEAN SURFACE FOR FIVE MODEL ATMOSPHERES 


Atmosphere 

Atmospheric conditions 

Pressure 

scale 

height, 

km 

Surface 

pressure, 

mb 

Surface 

temperature, 

°K 

Composition by 
mass, 
percent 

Composition by 
volume, 
percent 

Molecular 

mass, 

kg/kg-mole 

r, 

km" 1 

co 2 

n 2 

Ar 

C0 2 

n 2 

Ar 

Model I 

4 

262 

100 

0 

0 

100 

0 

0 

44 

-6.42 X 10“ 3 

13.38 

Model II 

4 

160 

100 

0 

0 

100 

0 

0 

44 

-4.59 X 10" 3 

8.11 

Model IH 

9 

219 

74.4 

12.8 

12.8 

68.5 

18.5 

13.0 

40.5 

-3.99 X 10" 3 

12.17 

Model IV 

20 

160 

25 

50 

25 

19 

60 

21 

33.6 

-3.08 X 10" 3 

10.77 

Model V 

20 

262 

| 25 

50 

25 

19 

60 

21 

33.6 

-5.23 X 10" 3 

17.56 
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Shadow of planet and atmosphere. Structure 
determined by angular width of Sun, and 
refraction, and Rayleigh scattering in planetary 



Figure 1.- Configuration for occupation. 
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Figure 2.- Configuration of source, planet, and observer. 
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Figure 3.- Differential refraction configuration. 
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(e) Mode! V; D = 8700 km. 
Figure 4.- Concluded. 
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(a) Mode! I; D = 8700 km. 
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(b) Model II; D = 8700 km. 

Figure 7.- The effect of Rayleigh scattering. 
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Figure 13.- Concluded. 
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